Collective effects in the dynamics of driven atoms in a high-Q resonator 
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We study the quantum dynamics of N coherently driven two-level atoms coupled to an optical 
resonator. In the strong coupling regime the cavity field generated by atomic scattering interferes 
destructively with the pump on the atoms. This suppresses atomic excitation and even for strong 
driving fields prevents atomic saturation, while the stationary intracavity field amplitude is almost 
independent of the atom number. The magnitude of the interference effect depends on the detuning 
between laser and cavity field and on the relative atomic positions and is strongest for a wavelength 
spaced lattice of atoms placed at the antinodes of the cavity mode. In this case three dimensional 
intensity minima are created in the vicinity of each atom. In this regime spontaneous emission 
is suppressed and the dominant loss channel is cavity decay. Even for a cavity linewidth larger 
than the atomic natural width, one regains strong interference through the cooperative action of a 
sufficiently large number of atoms. These results give a new key to understand recent experiments 
on collective cavity cooling and may allow to implement fast tailored atom-atom interactions as well 
as nonperturbative particle detection with very small energy transfer. 
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I. INTRODUCTION 

Cavity quantum electrodynamics CQED using opti- 
cal resonators has experienced important experimental 
progress in recent years. Several experimental groups 
achieved remarkable milestones in the realization of well 
defined strongly coupled atom-field systems [l|, 0, H, 0, 
i, S 0, i, S M HH, [H, [H, Q [H [Hi. This has lead 
to numerous applications like single atom trapping by a 
single photon [if, [3, 0] , conditional quantum phase shifts 
of very weak fields [3| , deterministic sources of entangled 
photons [s*, ^ and a single atom thresholdless laser [lol . 
Direct observations of the field mode structure [Tllll2l.ll8| 
and of the mechanical effects of the cavity field on the 
atomic motion [2, [^, [l9| have been reported culminating 
in the recent demonstration of cavity induced cooling of 
single trapped atoms [l|. 

Renewed interest was devoted to large ensembles of 
atoms commonly coupled to a cavity field with several 
modes, where collective atomic effects play a central role 
in the coupled atom-field dynamics [ll,[i3j[iM[3- These 
experiments have lead to unexpected results and opened 
new theoretical questions. As a particular example, col- 
lective phenomena in the presence of an external trans- 
verse driving field involving many atoms at different po- 
sitions are not fully understood [H, [l^ . 

In this work we study theoretically the dynamics of co- 
herently driven atoms in a resonator. Our investigation 
takes into account the atomic spontaneous emission, the 
finite transmittivity of the cavity mirrors and the spatial 
structure of the cavity mode. The system state is char- 
acterized by the atomic fluorescence rate and the signal 
at the cavity output and studied as a function of the 
system's parameter. To get further information on the 
system, we calculate the probe absorption spectrum and 
the field distribution in the vicinity of the atoms. 

Our results extend the studies of [l3| . There we showed 
that in the strong coupling regime the system dynamics 



exhibit enhanced cavity emission accompanied by sup- 
pression of fiuorescence which occurs, for more than one 
atom, when the atoms are spatially localized such that 
they emit in phase into the cavity mode. This phe- 
nomenon shares several analogies with the behaviour 
found in the case of a single atom inside a lossless res- 
onator [2^ [2l[. In fact, this behaviour can be traced 
back to destructive interference between the laser and 
the cavity field generated by atomic scattering, such that 
the atoms couple to a vanishing electric field. As a con- 
sequence, we show that the stationary cavity field is in- 
dependent of the number of atoms and cavity decay be- 
comes the dominant channel of dissipation. In a good 
cavity this allows to measure the light dissipated through 
cavity decay without destroying the interference, which 
is vital if one wants to get information on the cavity 
field and hence the current atomic positions. In addi- 
tion, when the strong coupling regime is achieved by a 
large number of ordered atoms the dynamics we find are 
consistent with the experimental observations by Vuletic 
and coworkers [l^, [lB| ■ 



This article is organized as follows. In section II the 
model is introduced. In section III we review the results 
obtained for a lossless resonator ^2fl|] and investigate the 
system's dynamics when the decay rate of the cavity is 
finite. In section IV the field inside the cavity is inves- 
tigated by means of (i) an additional weak laser coupled 
to the atom and (ii) an additional atom weakly coupling 
to the cavity field. In section V the scaling of the sys- 
tem dynamics with the number of atoms is studied and 
the results are discussed in connection with the experi- 
mental observations in [l^j [3] ■ In section VI the results 
arc summarized and discussed, and several outlooks are 
provided. The appendices report the details of the cal- 
culations presented in Sees. Ill and IV. 
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II. THE MODEL 




FIG. 1: N atoms couple to a ID optical resonator and are 
driven transversally by a laser which illuminate them homo- 
geneously. The inset shows the atomic transition which is 
relevant to the dynamics. 

We consider N identical and point-like atoms, whose 
dipole transitions couple resonantly with the standing- 
wave mode of an optical resonator. The atoms are as- 
sumed to be located along the axis of the resonator, which 
we denote with the x-axis, and their center-of-mass mo- 
tion is neglected. The relevant atomic degrees of freedom 
are the ground and excited electronic states \g), |e) of the 
dipole transition, which is at frequency loq- The transi- 
tion couples to the resonator's mode at frequency lOq and 
wave vector fc, and it is driven by a laser at frequency 
and Rabi frequency as shown in Fig. [TJ The laser 
is assumed to be a classical field. The dynamics of the 
composite system is described by the master equation for 
the density matrix p of atoms and cavity mode 

§-^P^^[H,p]+Cp + ^p (1) 

where H is the Hamiltonian for the coherent dynamics, 
and £, /C are the superoperators describing dissipation 
due to spontaneous decay and cavity losses. In the refer- 
ence frame rotating at the laser frequency wl the Hamil- 
tonian H has the form 

N 

H ^ -ndca^a~nAj2\(^)n{e\ (2) 

n=l 

N 
n=l 

where a, are the annihilation and creation operators 
of a cavity photon; (t„ = \g)n{e\, aj^ = |e)„((7|, are the 
dipole operators for the atom at the position a;„; Sc = 
LUi^ — LUc, A = LUi^ — luq are the detunings of the laser from 
the frequency of the cavity and of the dipole, respectively. 
The coupling constant between the dipole at position x„ 
and the cavity mode is g{xn) — gocoskxn, while the 
coupling with the driving laser depends on the atomic 
position through the phase 0„ = kxn cos 6, where 9 is 
the angle between the cavity axis and the propagation 
direction of the laser. Finally, the incoherent dynamics 



is described by the superoperators 

= i'^^^P^n - C^Icr«P - PCr^CTn) (3) 

n 

]Cp ~ —{2apa' — a'ap — pa'a) (4) 

where 7 is the rate of spontaneous emission of the dipole 
into the modes that are external to the cavity and k is 
the cavity decay rate. Collective effects in the sponta- 
neous decay are neglected here, as the average distance 
between the atoms is assumed to be of the order of several 
wavelengths. 



III. ENHANCED CAVITY EMISSION IN 
HIGH-FINESSE CAVITIES 

In [2^ it has been shown that the steady state of a 
lossless cavity, coupled to a dipole and driven transver- 
sally, is a pure state, such that the energy of the atom 
and cavity mode is conserved. This regime is accessed 
when the driving laser is resonant with the cavity mode. 
Then, the atom is in the ground state and the cavity 
mode field, generated by atomic scattering, is described 
by a coherent state whose amplitude is determined by the 
intensity of the laser. This can be easily seen in Eq. ([T]) 
for K = and iV = 1, after moving to the reference frame 
described by the unitary transformation 

V{(3) = exp(/3at - f3*a), (5) 

which corresponds to displacing the field inside the cavity 
by the amplitude /3 = -n&'^/g (here, g = g{x) ^ 22]). 
In this reference frame and for (5c = 0, Eq. ([2]) takes 
the form of the Jaynes-Cummings Hamiltonian and the 
steady state of the transformed master equation ^ is 
evidently the state \g,0)- In the original reference frame 
this corresponds to the steady state 

p,, = (6) 

which is a pure state. In fact, the state |g, (3) is eigenstate 
of H and is stable for k = 0. In particular, although the 
atom is driven both by laser and cavity mode, it is in 
the ground state as a result of the destructive interfer- 
ence between the atomic excitations induced by the two 
fields. Therefore, there is no atomic fluorescence and 
consequentely the only dissipation channel of the system 
is closed. State ^ exhibits thus the characteristic of a 
dark state. Moreover, as the atom does not scatter cav- 
ity photons and the cavity is assumed to be lossless, the 
energy of the cavity field is conserved. 

In any realistic setup an optical resonator has a finite 
decay rate n. When the atom is driven and 5c = 0, the 
field decay induces dephasing at the atomic position and 
the state ^ has a finite lifetime. However, it is reason- 
able to expect Eq. ^ to approximate the steady state 
for K sufficiently small. In the following, we investigate 
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the intensity of the fluorescence signal and of the sig- 
nal at the cavity output in different parameter regimes, 
thereby verifying under which conditions energy dissipa- 
tion through spontaneous decay can be neglected and 
when the state inside the cavity can be approximated by 
a coherent state. 

In this section we restrict to the case of one atom and 
take g{x) = g 0. The density matrix at time t for k ^ 
can be analytically determined by means of a perturba- 
tive expansion in the small parameter k, where cavity 
decay is assumed to be slower than the rate at which the 
atom reaches the steady state At this purpose, we 
rewrite the master equation ([T]) in the form 



d 1 / 



J p + k/CoP 



where 



hg{aa'^ + aV) + hVt{a^ + a) 



(7) 



(8) 



(A + i|)|e)(e| 



in— a' a 
2 



and Jp = japa\ JCqp — apa^ are the jump operators. 
The formal solution of Eq. ((7)) is [23| 



fulfilled for a wider range of values of k, as it is visible 
by comparing Fig[2^a) with Fig[2jb). In particular, for 
g = IO7 the signal /cav largely exceeds the fluorescence 
signal /at even for k > 7. 
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p{t) ^ Sit)piO) + dTS{t-T)iJ + KlCo)piT) (9) 

Jo 

with 5(i)/9(0) = exp{-iH^gt/n)p{0)exp{iHlgt/h). The 
perturbative expansion of Eq. ^ at second order in n 
is reported in Appendix A. The photon scattering rate 
at time t by the atom into the modes of the continuum 
is /at — 7Tr{CT^CT/9(t)}, and grows quadratically with k. 
At the cavity ouput the rate of photon scattering /cav = 
KTi{a^ap{t)} at time t and in lowest order is linear in k 
and /cav 2> /at- An instructive case is found in the limit 
A = and 3^7. Here, these expressions acquire the 
simple form 
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^cav ~ ' ^ 
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(10) 
(11) 



where Qi — 2g^ j^K is the cooperativity parameter per 
atom |25| . Thus, the two signals depend on n through the 
cooperativity parameter C\ and the factor kQF' /g^ , which 
is the decay rate of a cavity with mean photon number 
(n) = f^^/.g^. From Eqs. and ^ it is visible that 
for Ci ^ 1 the fluorescence signal is orders of magnitude 
smaller than the intensity at the cavity output. 

Figure [2] displays /at and /cav as a function of k and 
for two different values of g. The curves have been calcu- 
lated by solving numerically ([7|). Here, for a wide range 
of values of the cavity decay rate /cav exhibits a linear be- 
haviour as a function of k, while /at is quadratic. More- 
over, when g is increased (and thus when the coopera- 
tivity parameter is increased) the relation /at ^ /cav is 



FIG. 2; /cav (dashed line) and /at (solid line) as a function of 
K in units of 7. Here, 57 = 7, A = (5c = and (a) g{x) — 7, 
(b) g{x) = IO7. The orizontal dashed-dotted line gives the 
rate of fluorescence of the atom in free space. 

The zero-time correlation function g^^^ (0) of the signal 
at the cavity mirror gives further insight into the dy- 
namics of the cavity field. Figure ^h) displays g''^-' (0) 
as a function of k and Fig. [S^a) displays the correspond- 
ing average number of cavity photons (a^a). From these 
figures one sees that the cavity field exhibits a Poisso- 
nian behaviour for a fairly wide range of values of k, 
corresponding to large cooperativity parameters. This 
behaviour is verified even when the average number of 
cavity photons is very small (solid line in Fig. [3)Ja) and 
(b)). It shows that the cavity mode is in a coherent 
state, independently of the average energy of the cavity 
field. This behaviour contrasts dramatically with the an- 
tibunching observed when the pump is set directly on the 
cavity [11 HI]. 



IV. PROBING THE SYSTEM 

In this section we investigate the response of the system 
to a weak probe in the parameter regime for which the 
steady state is given to good approximation by state ©. 
We restrict to the case of one atom, whose dipole transi- 
tion couples to the driving field and to the cavity mode, 
and consider the spatial dependence of the coupling. 
The resonator's mode function is a standing wave with 
g{x) = gocoskx and the atom is assumed to be at po- 
sition X such that g{x) = g 0. The laser is a plane 
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B. The excitation spectrum is given by 
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FIG. 3: (a) Mean number of cavity photon and (b) second- 
order correlation function ^'•^'(O) as a function of k in units 
of 7. Here, A — Sc — 0, Q, = j and g{x) = 7 (dashed line), 
g{x) = IO7 (solid line). 



wave, and its phase 0„ depends on the atomic position 
Xn through the relation (j){x) — kxcosO, where 9 is the 
angle between the direction of propagation of the laser 
and the cavity axis. We assume that laser and cavity are 
resonant, and analyze the system's response to two types 
of probe: (i) a weak laser field, coupling to the atomic 
dipole, as shown in Fig. [H (ii) a second atom of a differ- 
ent species, whose dipole transition frequency is far-off 
resonance from the cavity frequency, thereby negligibly 
perturbing the system. 



A. Excitation spectrum 

We consider a probe driving the atom as illustrated in 
Fig. [31^ a) and evaluate the excitation spectrum, namely 
the rate of photon scattering into the modes external to 
the cavity, as a function of the detuning dp = up — ujl of 
the probe frequency up from the pump. For 5c = Q and 
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FIG. 4: (a) Excitation spectrum; A weak probe is coupled 
to the atomic dipole and its frequency is scanned through 
atomic resonance. The fluorescence signal is measured as a 
function of the probe detuning, (b) A second atom of another 
species, weakly coupled to cavity and pump fields, probes the 
cavity field inducing a position-dependent phase shift on the 
cavity field. The phase shift can be measured by means of 
homodyne detection or by measuring the fluorescence of the 
second atom. 

K = 0, the scattering rate of probe photons is evaluated 
for a probe Rabi frequency Jlp such that fip ^ g, 11,7. 
The details of the calculation are reported in Appendix 



w[5p) = 7?iJ7p 
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[5p[5p + /^)~-g^f + 5W/A 



(12) 



and it is plotted in Fig. [5] for two different values of the 
detuning A between atom and laser. From Eq. (fT^ it is 
evident that 'w{5p) vanishes at 5p = 0. This behaviour 
gives rise to a Fano-like profile of the excitation spec- 
trum as a function ol 5p [27i] . which is visible in Fig. [5l 
This profile is a manifestation of destructive interference 
between the excitation paths contributing to the atomic 
dynamics, which can be identified with the absorption 
of photons from the laser and from the cavity field [28[. 
Interference is at the origin of the two resonances visible 
in Fig [5] where the rate of photon scattering is maxi- 
mum. They correspond to values of the probe detuning 
5p ^ 5±, with 



and have width 7-1- , which for 
simple form 




g2 > take the 



(13) 



In the strong coupling regime these resonances corre- 
spond to the dressed states of the atom-cavity system, 
and their widths determine the characteristic time-scales 
of the system's dynamics. Thus, for k 7^ enhanced cav- 
ity emission accompanied by suppression of fiuorescence 
are achieved when min(7+,7_) > n. 

It is remarkable that w{5p) does not depend on Vl and 
thus does not depend on the average number of photons 
inside the cavity. In particular, the position and width 
of the resonances are the ones found for an atom in an 
empty cavity. This result can be simply explained by ob- 
serving that the field at the atomic dipole vanishes. In 
other words, in the reference frame described by the uni- 
tary transformation ([5]) the absorption of a probe photon 
induces a transition |(7,0) — > |e, 1), whereby |e,l) is the 
superposition of the eigenstates of the system at frequen- 
cies 5± corresponding to the resonances of the excitation 
spectrum. For k ^ we have verified that also the split- 
ting at the cavity output depends on g{x) and not on SI. 
This property suggests a use of this interference effect 
in order to probe the atomic position inside the cavity 
without significantly perturbing the system. 



B. A second atom probing the cavity field 

The electric field inside the cavity can be probed by 
means of an atom which is weakly coupled, as shown in 
Fig. IHb). This can be, for instance, an atom of other 
species whose dipole transition frequency is far-off reso- 
nance from the cavity and the driving laser. This atom 
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FIG. 5: Excitation spectrum w{5p) as a function of the probe 
detuning Sp in units of 7, for the parameters g{x) = y, 5c = 0, 
K = and (a) A = 0, (b) A = -27. 



thus experiences a small a. c. -Stark shift (5atom, whose in- 
tensity is a function of the distance d = x' — x from the 
atom which pumps the cavity. It can be measured with 
an homodyne detection of the cavity output field, where 
the pump field is the local oscillator, or by measuring the 
fluorescence of the probe atom. The shift iJatom is plotted 
in Fig. [6] as a function of the distance d for different val- 
ues K. We observe that cavity decay tends to cancel the 
spatial modulation of the total electric field, which never 
vanishes inside the cavity for k ^ 0. If the mechanical 
effects of light are considered, then Fig. [6] corresponds to 
the potential that the probing atom experiences. Hence, 
the latter may feel a binding or repulsive force in the 
vicinity of the first atom, depending on the sign of the 
detuning A2 of the laser from the probing atom reso- 
nance. 
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FIG. 6: a. c. -Stark shift (5atom of the probing atom as a func- 
tion of the distance from the atom pumping the cavity, which 
is at an antinode of the standing wave. Here, go = = 7, 
A = (5c = 0, and the laser propagation direction is perpen- 
dicular to the cavity axis. The detuning of the probe atom 
from the cavity frequency is A2 = IOOO7. Solid line: k = 0; 
Dashed line: k = 7; Dash-Dotted line: k = 27. 



V. SCALING WITH THE NUMBER OF ATOMS 

A. Two atoms inside the resonator 

So far we have considered that only one atom couples 
to the cavity mode. If a second atom of the same kind 
is inside the cavity and is illuminated by the laser, the 



dynamics is in general non-trivial. For k = suppression 
of fluorescence is observed when the atoms are at a dis- 
tance Ax which is an integer multiple of the wave length 
A. When this occurs, the two atoms are coupled with 
the same coupling constant g to the cavity mode, and 
from Eq. ([T]) it can be verifled that the state jgi, 32, /?) is 
the steady state of the system with (3 = —^/g. Hence, 
the total electric field vanishes at both atoms, while the 
cavity field is the same as when only one atom couples 
to the resonator. In general, one can define the function 



/3(x) = n{x)exp{i{TT + kxcos9))/g, 



(14) 



where g — g{x) ^ and which includes inhomogeneity of 
the pumping field. Then, the condition for suppression 
of fluorescence with two atoms is fulfilled whenever two 
positions exist such that f3{xi) = P{x2), where the atoms 
are located. Clearly, there may exist parameters regimes 
for which function (jl4p is not periodic and a non-trivial 
solution for suppression of fiuorescence with more than 
one atom does not exist. 

Figures [7] and [8] display the average number of photons 
and the excited state populations as a function of the 
relative distance between the atoms, assuming that one 
atom is fixed at the antinode of the standing wave and 
that both atoms are homogeneously driven by the laser, 
which propagates perpendicularly to the cavity axis, i.e. 
6 = tt/2. Here it is visible that the excited state popula- 
tion of both atoms vanishes at X2 = 0, A. At these points 
the field inside the cavity is different from zero, and it 
is a local minimum as a function of X2, as it is particu- 
larly evident in Fig [71 The population of the first atom 
vanishes as well when the second atom is at a node of 
the standing wave. In this case, the dynamics of the cav- 
ity are determined by the coupling with the atom at xi , 
whereas the population of the second atom is determined 
by the laser intensity, as if it were in free space. 
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FIG. 7: (a) mean number of photons and (b) excited state 
population of atom at xi (solid line) and at X2 (dashed line) 
as a function of X2- Here, xi is an antinode of the standing 
wave g{x) — go cos kx. The parameters are k = O.27, S7 = 7, 
go — IO7, A = IOO7, Sc = 0, and the laser propagation 
direction is perpendicular to the cavity axis. 

Note that Fig. [7] displays the situation when the atoms 
are driven below saturation. Here both scatter coher- 
ently into the cavity mode and a second type of interfer- 
ence effect occurs for X2 = A/2: At this point the cavity 
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field vanishes, wiiile the atomic populations are equal and 
different from zero. In fact, for X2 = A/2 and below satu- 
ration the atoms scatter coherently into the cavity mode 
with opposite phase. At saturation, on the other hand, 
the scattered light is mostly incoherent, and the cavity 
field does not vanish at this point, as shown in Fig. [51 
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FIG. 8: Same as Fig.[7l Here, k = O.OI7, = 7, 5 = IQ7, A = 
5c ~ 0, and the laser propagation direction is perpendicular 
to the cavity axis. 

In Fig. [9] the ratio between the cavity and fluorescence 
signal is displayed for the case illustrated in Fig. [7| when 
the atoms are driven below saturation and the laser is 
orthogonal to the cavity axis. Here one sees clearly that 
this ratio is maximum when the atoms are a wavelength 
apart, namely where the function (|14|) assumes the same 
value at the atomic positions. For k ^ absolute maxima 
are found when the atoms are at the antinodes of the 
cavity mode, where the cooperativity parameter is largest 
and the total electric field vanishes. 
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B. N atoms inside the resonator 

The dynamics of the coupled system for generic pa- 
rameters and number of atoms are very complex. Nev- 
ertheless, insight can be gained in the limit in which the 
atoms are driven below saturation. This assumption en- 
ables one to adiabatically eliminate the atomic degrees 
of freedom from the cavity equation, and corresponds to 
the parameter regime I7/2 -t- iA| 3> VNg, ^/Nfl, when 
the collective dipole is driven below saturation. Under 
these conditions the density matrix of the field pt obeys 
the equation 



dt 2^' 



k) {2apfa' — a' apt — pfO^ a} (15) 
+iS'[a^a,pt]-i[i^a^ +Ca),pi] 



where j'{N) — Nsj is the cavity decay rate due to 
photon scattering by spontaneous emission, S'{N) = 
6c — NsA contains the a. c. -Stark shift due to the medium, 
and ^ is the cavity drive mediated by the dipoles. 



.7\ E„s(^n)f^e^'*" 



(16) 



Here, s = ^n/N where s„ is defined for the atom n as 



(7/2)2 + A2 



(17) 



From Eq. p^)) it is visible that the system dissipates with 
rate 7' + ^, which determines the rate at which the steady 
state is reached. The steady state of (fTS)) is pf^gg ~ \a){a\, 
where \a) is a coherent state with amplitude 



(7'(iV) + K)/2 + i(5'(iV) 
E„ff(x„)e'*" (7/2 + iA) 



E„l5(a;„)P {Ns^ + K)/2-i{5c-NsA) 

(18) 



FIG. 9: (a) Ratio F = /cav//at for two atoms as a function 
of their position x\ and X2 inside the cavity for the same 
parameters as in Fig. [T] (b) Ratio as in (a) as a function 
of X2 for xi — 0. The dashed line shows F for the same 
parameters but k — 'y. Note that xi, X2 are plotted modulus 
A, and x\ ^ X2- 

These considerations can be extended to three dimen- 
sions in a straightforward way. The three-dimensional 
pattern is found taking into account the phase of the 
pump, which we always assume to be orthogonal to the 
cavity axis. The zeros of the electric field are then dis- 
tributed according to a Body-Centered-Cubic lattice with 
distance A/2 between adiacent planes [H, |2^. Fluores- 
cence is suppressed when the atoms are localized at these 
points, thus forming a stationary pattern. 



and which is the sum of the electric fields scattered at 
each atom. In fact, in this regime the collective dipole is 
driven well-below saturation and radiation is scattered 
elastically into the cavity mode. 

For a large number of atoms the contributions of each 
atom sum up so that the field amplitude a exhibits a 
narrow peak at the maximum value ao as a function 
of the mean square deviation of the phase A0„. This 
maximum corresponds to the case when all atoms scat- 
ter in phase, namely when they are distributed at the 
points {xi, . . . , xjv} where the function (jl4p acquires the 
same value. The necessary condition that this situation 
is verified is that j3(x) is periodic, as we have previously 
observed. In the following we assume that the laser prop- 
agates perpendicularly to the cavity axis, i.e. 9 = 7r/2. 
Hence, I3{x) has periodicity equal to A. Assuming that 
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the atoms scatter in phase into the cavity mode, the am- 
phtude of the cavity field is given by 

n Ns{j/2 + iA) 
° 5 7Vs(7/2 + iA) + K/2-i(5c ^ ' 

where g = g(xi) = ... = gix^) and g 7^ 0, namely 
the atoms are spatially distributed in a pattern which 
has spatial periodicity A. The pattern is here assumed to 
have low filling factor, such that sub- and supperradiance 
effects in the scattering in free space are negligible. Note 
that for large filling factors subradiance may give rise to 
other meta-stable states of the collective dynamics. 

1. Stability of the atomic patterns 



in first order in 6xn, nor does the force for small fluctu- 
ations in a. Remarkably, Eq. (|2ip is independent of A. 
Moreover, its intensity depends on the mean number of 
cavity photons, and thus on the pump intensity. This 
result has been confirmed by numerical simulations and 
is in line with the experimentally observed dependence of 
enhanced cavity emission on the intensity of the pump, 
showing that the effect manifests itself when the pump 
intensity exceeds a threshold value [T^ . We remark that 
result (|2ip is valid when the ratio 6c/N is sufficiently 
small, so that to good approximation the field inside the 
cavity is given by = {—lY'^^fl/ go and the total field at 
the atomic positions almost vanishes. The dependence of 
the system dynamics on the atom number N is discussed 
in the following subsection. 



The atomic pattern in a standing wave cavity is invari- 
ant per translation by A. In principle, there is an infinite 
number of patterns for any value of x in the range [0, A), 
corresponding to different values of the coupling constant 
g{x). When the mechanical effects of light are taken into 
account, however, the equilibrium positions are at the 
antinodes of the cavity-mode standing wave a; = or 
X — X/2. At these points, in fact, the force on the atoms 
vanishes. This is evident when we consider the force Fn 
entering the semiclassical equation of the motion for the 
atom at Xn moving along the cavity axis 

Fn = hkUo\a\'^ sm{2kxn) 

+2hklm{r]*ffa} sin(fca;„) (20) 

where Uq = gQA/{A'^ + 7^/4) is the light shift due to 
the coupling to the cavity. To = gl^/2{A'^ + l'^ 1^) the 
rate of dissipation, and 77off = ^go / {—lA + ^ / 2) the term 
due to the transversal pump on the atoms. The param- 
eter a describes the field amplitude, which evolves ac- 
cording to Eq. IjlSp. We denote with even (odd) pat- 
tern the atomic pattern where the atoms are localized 
at the equilibrium points Xn^ such that cos 
(cos kxn'' = — I). For k = and 5c = 0, one has suppres- 
sion of fluorescence when the cavity field coherent state 
has amplitude a = = f7/goe'(^+^^'', with £ = 0, 1 de- 
pending on whether the pattern is even or odd. Hence, 
the cavity fields due to each pattern differ by a phase tt. 
Numerical studies have reported selforganization of the 
atoms in these patterns [29|. Stability is found for small 
(but non- vanishing) negative values of 5c, as we have ver- 
ified numerically . From Eq. ([^0]) we can estimate the 
force 5fn around these points when the atoms undergoes 
a small displacement 5xn from the equilibrium position 
Xn^ . At first order in 5xn the force takes the form 

5U^2ne(-\ ^fSxr. (21) 

and it is clearly a restoring force for 5c < This condi- 
tion is sufficient, since the field amplitude a does not vary 



2. The cavity field when the atoms emit in phase 

We now assume that the atoms are localized in an 
even pattern, such that the cavity field amplitude is 
Pi = /?o = —^/goi and disregard the mechanical effects 
of light. For k — Q and 5c — we recover from Eq. ([TO]) 
the result ao = /3o- Thus, in this limit the stationary field 
is independent of the number N of atoms and of the de- 
tuning A between laser and dipole transition. The field 
amplitude achieves the maximum value as a function of 
5c for 5' = 0, corresponding to the condition 5c = NsA. 
This is visible in Fig. I10[ where the average number of 
photons is plotted as a function of 5c and for two different 
values of A. For A Q the detuning 5c = NsA is the 
a. c. -Stark shift of the cavity mode frequency due to the 
coupling with the atomic dipoles. For this value the clas- 
sical field drives the system resonantly, and the amount 
of energy transferred into the cavity mode is maximum. 
Note that the position of the resonance 5c — NsA scales 
linearly with N and for |A| ^ 7 is inversely proportional 
to A. The corresponding linewidth 'y' + k scales as N/ A^ 
for 7' » K. Obviously, large values of k broaden the res- 
onances. 
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FIG. 10: Mean number of cavity photons as a function of 
5c in units of 7. The solid line corresponds to k = 0, the 
dotted line to k = O.OI7. For N = 1 atom, the parameters 
are n = 5 = O.I7, and (a) A = 0, (b) A = 7. 

At the amplitude of the cavity field of Eq. (fTO|) the 
population of the excited state of an atom in any of the 
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FIG. 11: Excited state population as a function of 5c in units 
of 7, Same parameters as in Fig. 1101 



pattern positions is given by 



(7/2)2 + A2 (y + k)2/4 + (5'2 



(22) 



and is displayed in Fig. [TT] as a function of Sc for some 
parameter regimes. Clearly, for k = and 6c = the ex- 
cited state population vanishes, indicating that the atoms 
stop fluorescing. For A the population lie exhibits 
a maximum, which is located at Sc NsA{l + j'^/AA'^) 
for |A| ^ 7. For k ^ the center-frequencies of the 
resonances are shifted by an amount proportional to the 
cavity decay rate, the curves are broadened, and the ex- 
cited states population does not vanish at Sc = 0. 

The behaviour of the system as the number of atoms 
N is varied exhibits remarkable features. In fact, N ap- 
pears in the denominator of Eqs. and (P^ . scaling 
the atomic effects in the cavity dynamics. In particu- 
lar, a critical value No for the number of atoms can be 
identified, such that for N ^ Nq the coupling with the 
atoms affects relevantly the cavity dynamics, whereas for 
N Nq atoms and cavity are weakly coupled. 

For A — Sc ^ one finds the value A'o = k/s7 = 1 /2Ci 
where Ci = 2(7g/K7 is the one- atom cooperativity param- 
eter [25|. Thus for N Nq the system is characterized 
by a large cooperativity parameter. In particular, when 
N No the excited states population in Eq. ([^ ac- 
quires approximately the value as in free space, while 
the cavity field amplitude scales linearly with the num- 
ber of atoms. There is thus no back-action of the cav- 
ity on the atomic dynamics, since the cavity decay rate 
is faster than the rate at which the atomic degrees of 
freedom reach their steady state. On the other hand, 
when N ^ No the field amplitude tends to the asymp- 
totic value a — > —il/g, while HelgOc /N'^. Thus, the 
power dissipated by spontaneous emission scales with 
1/iV, while the signal at the cavity output is constant and 
independent on the number of atoms. Figure [T^] displays 
the signal at the cavity output and the total fluorescence 
signal evaluated from Eqs. and ([^ . respectively as 
a function for TV. For these parameters No ^ 10"^. 

An analogous behaviour can be found for large val- 
ues of A, and is illustrated in Fig. [T31 For |A| 3> 7, k 
the critical number of atoms, determining the regime of 
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FIG. 12: 7cav and total 7at as a function of A'^ for Q, 
K = 10"^7, and A = 5^ = 0. 
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FIG. 13: Icav and total /at as a function of A'^ for Q. = g 
K. = IO7, A = -IOOO7 and (a) 5c = 0, (b) 5c = -67 



strong coupling, is given by A'q.a = l^l'^/ff^i and en- 
hanced cavity emission accompanied by suppression of 
fluorescence is observed for N ^ iVo,A- This behaviour 
is also found for values of the detuning Sc 7^ 0, as shown 
in Fig. [T3Fb) for Sc ~ —87, provided that N is suffi- 
ciently large to fulfill the relation \Sc\ ^ iVs|A|. Note 
that for values of N such that Sc = iVsA, namely when 
a collective state of the system is driven resonantly, the 
signals in Fig. fTSFb) exhibit a maximum. Nevertheless, 
as N increases, /cav tends asymptotically to the value 
lew K|il/goP which is independent, among others, of 
N and of A. It should be noted that in this case increas- 
ing N corresponds to increasing the width of the window 
around the value (5c = appearing in the atomic pop- 
ulation as a function of Sc^ as shown in Fig. [TlJ Thus, 
the condition |(5c| <C iVs |A| corresponds to values of the 
detuning Sc which are much smaller than the a. c. -Stark 
shift, hence for which the condition of destructive interef- 
erence is still (although approximately) fulfilled. 

The parameter regimes discussed in Fig. I13f b) are con- 
sistent with the ones of the experiment by that 
reported a rate of emission into the cavity modes ex- 
ceeding by orders of magnitude the rate of fiuorescence 
into the modes external to the cavity. This observa- 
tion was accompanied by the measurement of a coherent 
cavity field whose characteristic gave evidence of atomic 
self-organization. This behaviour has been explained as 
Bragg scattering of the pump light by the atomic grat- 
ing. However, from the results presented in this section 
we can argue that Bragg scattering is actually suppressed 
in this regime, as the experimental regime of hj can be 
classified to be in the region with N ^ No- In fact. 
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in the strong couphng regime a cavity field establishes 
when the atoms organize spatially, that cancels out with 
the pump at the atomic positions. As a consequence the 
atoms decouple from the cavity and pump field, and are 
in the ground state. Therefore, there is no fluorescence 
nor superradiant scattering into the cavity mode. In this 
regime the main source of dissipation is through cavity 
decay. We remark that these dynamics is encountered 
also in the case of a single atom. In fact for strong cou- 
pling the stationary cavity field is solely determined by 
pump intensity and cavity coupling at the atomic posi- 
tion, while the atoms are in the ground state. In par- 
ticular, the regime of large N in Figs . fl2lfT3l corresponds 
to large cooperativity parameters, while in Fig. [2] strong 
coupling is achieved for small k. Note that, differently 
from optical bistability [33] , where bistable dynamics are 
observed for large cooperativity, here there is only one 
steady state, where the atoms are in the ground state. 

In summary, the enhanced cavity emission of [l6j can 
be traced back to an interference effect between pump 
and cavity field, which is established for large coopera- 
tivity parameters. Although in this treatment we have 
neglected the center-of-mass motion, this hypothesis is 
supported by the stability of the pattern in this param- 
eter regime. 

Finally, it should be observed that one important con- 
dition for these dynamics is that the atoms are localized 
according to a stationary pattern. This condition consti- 
tutes a substantial difference to the collective scattering 
via acceleration observed in the dynamics of the collec- 
tive atomic recoil laser [l^ US] . 



VI. DISCUSSION AND OUTLOOK 

Collective effects play an important role in the dynam- 
ics of N coherently driven atoms within a cavity mode. 
For a resonant laser and a sufficiently large cooperativ- 
ity parameter the atomic scattering of photons into the 
cavity field may exceed the scattering into the free-space 
modes by several orders of magnitude despite weak cou- 
pling and low excitation of each individual atom. In 
this regime the cavity output exhibits Poissonian photon 
statistics independent of the mean cavity photon num- 
ber. In addition the probe excitation spectra of the atoms 
reveal the cavity vacuum-Rabi splitting even for strong 
pump fields, when the mean number of cavity photons is 
large. 

The phenomenon can be understood as interference be- 
tween the transverse pump and the cavity field, which is 
established when the coupling between atoms and cavity 
mode is sufficiently stronger than other effects determin- 
ing the system's dynamics. Remarkably, the establishing 
of this regime corresponds to the situation in which the 
total electric field at the atomic position vanishes. For 
two or more atoms inside the cavity these conditions are 
accessed when the atoms are distributed according to a 
spatial pattern with periodicity equal to the mode wave- 



length. By means of a simple model we have shown that 
stable patterns are achieved for suitable laser and cavity 
parameters, when the locations of the pattern are at the 
antinodes of the cavity standing wave. We have identified 
two patterns, which correspond to fields inside the cavity 
which are shifted by a phase tt. The results predicted 
by this model are in qualitative agreement with the dy- 
namics reported in jl^, iXQ, i29| , and provide a physical 
picture of the phenomena observed. It should be noted 
that the cavity used in flEi [l6j is multimode, whereas in 
this work we consider a single mode cavity. Nevertheless 
the dynamics reported in [l6| can be reproduced with a 
model consisting in a single mode cavity and two level 
atoms, showing that the basic physical phenomena can 
be traced back to the interference effect discussed here. 

The phenomenon of interference in the driven Jaynes- 
Cummings model has been denoted with "cavity induced 
transparency" by Rice and Brecha ^2§\, and it can be 
traced back to the classical dynamics of two coupled 
damped oscillators [1^ . Rice and collaborators have the- 
oretically investigated the dynamics of a classical dipole 
coupling resonantly to a cavity mode when the cavity is 
driven p6l Issj . In this case the field due to atomic po- 
larization cancels out with the drive on the cavity. Due 
to this effect the cavity electric field vanishes. Thus, the 
phenomenon is established when the dipole decay rate is 
smaller than the cavity decay rate and, once this regime 
is accessed, energy is dissipated mainly by spontaneous 
decay. This situation might seem equivalent in many re- 
spects to the case discussed by Carmichael and coworkers 
in [13], where the role of cavity and atom are exchanged. 
Nevertheless, when the cavity is driven quantum noise 
and saturation effects on the dipole give rise to devia- 
tions from its classical behaviour and thus from inter- 
ference. Here interference is recovered for a sufficiently 
large number of dipoles N, so that the collective dipole is 
to good approximation an oscillator. Another interesting 
difference between the driven-cavity and the driven-atom 
case is the signal at the cavity output. When the cavity 
is driven and 7 <C k, the g^^-* (t) function is antibunched 
at r = pol . [35I ] . On the contrary, when the atoms are 
driven and k <C 7, we have shown that ^^^•'(O) = 1 even 
when the mean energy of the cavity field energy is very 
small. 

It is instructive to compare the phenomenon of sup- 
pression of fluorescence investigated in this work with 
the phenomenon of electromagnetically induced trans- 
parency manifesting itself in driven multilevel atomic 
transitions [s^]. The two types of interference arise be- 
cause of different dynamics: In EIT the atomic polar- 
ization is orthogonal to the field polarization, so that 
the atom does not absorb photons. In " cavity- induced 
transparency" the laser and the cavity field cancel out, 
so that the total electric field at the atom is zero. It is 
this very property that leads to the vacuum Rabi split- 
ting observed in the excitation spectrum even when the 
mean energy of the cavity field is significantly large. 

There are several interesting questions which are worth 
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investigating starting from these results. For instance, do 
other patterns exist, than the ones found, which may be 
stable and bring to a different steady state of the cav- 
ity field? In this context, we have considered a transverse 
pump whose propagation direction is orthogonal to the 
cavity axis. In this way, the periodicity of the pattern is 
determined solely by the periodicity of the cavity stand- 
ing wave. When the angle between cavity and pump is 
different the situation may change drastically, even giv- 
ing no a priori possibility of finding a stable pattern for 
more than one atom. For other cases, when the cav- 
ity mode is not a standing wave, but, say, a ring cav- 
ity [H, again other dynamics are expected. Such 
questions will be tackled by treating systematically the 
mechanical effects of light-atom interaction, and will be 
subject of following works 32]. 

For systems of one or few atoms in high-Q cavities 
the interference phenomenon presents several potentiali- 
ties for implementing coherent dynamics of quantum sys- 
tems. For instance, the vacuum Rabi splitting observed 
by probing the system depends on the position of the 
atom in the mode, and may allow to determine the spa- 
tial mode structure [lH , as well as to implement feedback 
schemes on the atomic motion [s^ |3^. Moreover, sev- 
eral experimental setups can presently trap single or few 
atoms and couple them in a controlled way to the cav- 
ity field d, [a, [m [II] . The dynamics discussed here 
can be applied for instance to implementations of quan- 
tum information processing, since interference effects are 
rather robust against noise and decoherence. In addi- 
tion, the coherence properties of the transmitted signal, 
which are preserved even for very small photon numbers, 
suggest an alternative kind of photon-emitters to the one 
investigated in [1, [1, [13] . 
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APPENDIX A: EVALUATION OF THE 
PERTURBATIVE CORRECTIONS IN n 

In this appendix we present the main steps to de- 
termine p{t) in ^ to second order in k. In the dis- 
placed frame defined by the unitary transformation ([S]) 
p[t) = V{(3yp{t)V{/3), and Eq. (0) takes the form 

p{t)^S{t)m+ I dTS{t^T){KlCo+-/J)p{T) 

Jo 



where p{0) = 2?(/?)tp,,2?(/?) = \g,0){g,0\ and 

o /n\^ 

ICoX = aXa^ (aX + Xa^) + { -] X 

9 \gj 
S{t)X = U,s{t)XU,s{ty 

Ucifit) = exp (^-^Hcst^ 

Here, ICq, S{t) are the transformed superoperators de- 
fined on a density matrix X and Ucs{t) is a non-unitary 
operator, with Hgg = T>{(3y Heg'D{f3). The operator HeS 
is non-Hermitian. Its spectrum and eigenvectors are cal- 
culated by solving the secular equations for the right and 
left eigenvectors of HcS, according to Hcs\v\) = Ti\\v\)^ 
hI^\v\) = h\*\v\). The states {\v\), \v\)} constitute a 
biorthogonal basis, such that 1i'a)(wa| = 1- In this 
basis the operator ?7off (t) can be written as 

C/eff(i)=Ee"'^>A)M (Al) 
A 

We expand now Ucs{t) in second order in the parame- 
ter K, U,s{t) = U^^{t) + nU^s\t) + n-'U^lHt), where the 
superscript indicates the corresponding order in the per- 
turbative expansion. In order to evaluate these terms we 
define 

Hcff = i?iff^ + kV (A2) 

with 

= ?ig(aa^ + aV)-?i(A + i|) le)(el (A3) 
V = --[a^a+-J+--ia^ + a) (A4) 

and solve the eigenvalue equation at second order in k, 
obtaining 

A = A(°)+«A(i)+k2a(2)+o(«;3) 

and analogously for the left eigenvectors, where JV = 
In particular, the solutions at zero order have 

the form 

Ai"L = (^A + i^ T y(A + i^)' + Vnj (A5) 

with the respective right eigenvectors 

\v'^l) = a„,± Je, n - 1) + bn,±\g, n) (A6) 

and 

bn,± = gVnAfn,±^ 
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while the left eigenvectors \v) have the form |w„.±) = 

a* _|_|e,n-l) + 6* whereA/'n,± = {vn,±\vn,±)- Us- 

ing this expansion, the terms of the perturbative expan- 
sion of the operator J7eff(t) are immediately found, and 
the evaluation of p{t) = p'^°\t) + Kp'^^\t) + n'^p'^'^^t) + 
o{k^): consists in the evaluation of the integrals 

p(«)(t)=5(°)W|5,0)(.9,0| (A7) 

{t) = /* dr5(") {t - r)/Co5(o) (r)|.g, 0) (g, 0| 



p(2)t)=5(^)(t) 1.9,0) (.9,01 



dr 



5(")(t~T)J5(2)(r) 



{t - r)/Co5(i) (r) + S^^^ {t - t)/Co5(") (t)] \g, 0) (5, 

where S^^\t)X ^ J2l=o U^cS ^\t)XU^p {t^ . Note that 
in (|A7|) we have omitted to write the terms containing 
JS'^^\T)\g, 0){g, 0| (with j = 0, 1), since they vanish. 



We assume that at t = the system is in the stationary 
state of the driven dipole and cavity system, the probe 
field is a coherent state of amplitude 77 such that jjyp <C 1 
and the other modes of the electromagnetic field are in 
the vacuum |0k)- As we are interested in the probability 
that a probe photon is scattered into the modes of the 
e.m.f.-field, the initial state is given with probability \ri\'^ 

by 

|^,) = |g,/3;lp,0k) (B4) 
and it is at energy Ei = hA + hSp, while the final state 



IV/.k) = \g,(3] Op, Ik) 



(B5) 



The transition amplitude is the element of the scattering 
matrix 5,./^, 

5,;/, = ~2in hm S^^^ [Ej, - E,)T,f, (E,) (B6) 



APPENDIX B: EVALUATION OF THE 
EXCITATION SPECTRUM 

We evaluate the excitation spectrum by calculating the 
transition amplitude which describes the scattering of a 
probe photon into the modes of the electromagnetic field, 
into which the dipole spontaneously emit. The modes of 
the electromagnetic field are here treated quantum me- 
chanically. The Hamiltonian determining the dynamics 
is 



H' — H + -ffpiobo 



(Bl) 



where H is defined in Eq. iJprobo = fiSpb''b + V, with 
6, 6^ annihilation and creation operators of a probe pho- 
ton, Sp detuning of the probe from the cavity frequency, 
and 



V = hilp {bai + b'<a„ 



(B2) 



the interaction of the probe with the dipole, with Qp 
vacuum Rabi frequency The term i?cmf describes the 
coupling of the dipole to the other external modes of the 
electromagnetic field , Hcmf — X^k ^'^k^k^k + W, where k 
labels the mode at frequency (5k (in the reference frame of 
the drive) and wave vector k, with corresponding creation 
and annihilation operators fej^., fek, and W describes the 
interaction with the atomic dipole. 



= ^ %k (cT^bi, + abl 



(B3) 



Here, is the vacuum Rabi frequency for the coupling 
of the mode to the dipole. 



where (5^^-' {E) is the diffraction function, 

1 sm(CT/2n) 

TT E 



(B7) 



and Tfi{Ei) is the transition matrix element, which at 
lowest non-vanishing order has the form 



Ei — Heff 



-V\^,) (B8) 



Here, H'^^ is the effective Hamiltonian, 



H's = H,s + hSpb^b 



(B9) 



where Hcs is given in ([8|). At lowest order in the 
transition matrix element (IB8I) has the form 



1 



T,^,(£;,) = ?i"5k^^p(e,/3;0p,0k|- pj^le,/?; Op, Ok) 



(BIO) 



where Clp — r/flp. The solution of (|B10[) can be easily 
found in the reference frame defined by ([5]): Here, H'^f^ — 

H^g^ + hSpb'fb, where we have used (IA3[) . Finally, we 
obtain 



(5p[5p + A + 17/21-52 



(Bll) 



By substituting this result in Eq. (jB6[l we find the tran- 
sition amphtude. The rate is found from (|B6p after 
summing over all modes of the continuum and taking the 
modulus squared divided by the time T [2^ . 
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